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ABSTRACT 

We present new fully self-consistent models of the formation and evolution of isolated dwarf 
galaxies. We have used the publicly available N-body/SPH code HYDRA, to which we have 
added a set of star formation criteria, and prescriptions for chemical enrichment (taking into 
account contributions from both SNIa and SNII), supernova feedback, and gas cooling. We 
extensively tested the soundness of these prescriptions and the numerical convergence of the 
models. The models follow the evolution of an initially homogeneous gas cloud collapsing 
in a pre-existing dark-matter halo. These simplified initial conditions are supported by the 
merger trees of isolated dwarf galaxies extracted from the milli-Millennium Simulation. The 
star-formation histories of the model galaxies exhibit burst-like behaviour. These bursts are a 
consequence of the blow-out and subsequent in-fall of gas. The amount of gas that leaves the 
galaxy for good is found to be small, in absolute numbers, ranging between 3 x 10 7 Mq and 
6 x 10 7 Mq. For the least massive models, however, this is over 80 per cent of their initial 
gas mass. The local fluctuations in gas density are strong enough to trigger star-bursts in the 
massive models, or to inhibit anything more than small residual star formation for the less 
massive models. Between these star-bursts there can be time intervals of several Gyrs. 

The models' surface brightness profiles are well fitted by Sersic profiles and the correla- 
tions between the models' Sersic parameters and luminosity agree with the observations. We 
have also compared model predictions for the half-light radius R c , central velocity disper- 
sion c c , broad band colour B — V, metallicity [Z/Zq] versus luminosity relations and for 
the location relative to the fundamental plane with the available data. The properties of the 
model dwarf galaxies agree quite well with those of observed dwarf galaxies. However, the 
properties of the most massive models deviate from those of observed galaxies. This most 
likely signals that galaxy mergers are starting to affect the galaxies' star-formation histories 
in this mass regime (M > 10 9 Mq). 

We found that a good way to assess the soundness of models is provided by the combi- 
nation of R c and er c . The demand that these are reproduced simultaneously places a stringent 
constraint on the spatial distribution of star formation and on the shape and extent of the dark 
matter halo relative to that of the stars. 
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1 INTRODUCTION 

Dwarf galaxies (DGs), galaxies with blue absolute magnitude 
— 14 > Mb > —18 m ag, are the most numerous t ype of galaxies 
in the nearby universe jFerguson & Binggel i 1994). Besides their 
number there are various reasons why DGs are very important in 
modern astronomy, most notably their cosmological relevance. 

In the Lambda cold dark matter (ACDM) cosmological model, 
present day galaxies are the product of a series of mergers. The 
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early universe consisted of a uniform density, superimposed with 
small Gaussian fluctuations. As their amplitude increases monoton- 
ically toward smaller scales, they give, after inflation, rise to a myr- 
iad of small virialized dark matter (DM) haloes. These small haloes 
collapse first under influence of their own gravity, followed by 
the collapse of haloes with larger masses, swallowing the smaller 
haloes. We thus have hierarchical clustering, where massive galax- 
ies, such as our own Milky Way, have formed through the merger 
of DGs. Dwarf galaxies are therefore an important link in the chain 
of hierarchical clustering. Moreover due to their relatively low total 
mass (0(1O 7 — 10 9 Mq)) they are quite sensitive to the effects of 
gravitational interactions, star formation (SF) and supernova (SN) 
explosions. On the one hand DGs are therefore the ideal objects to 
study when researching SF and SNe in a galaxy environment, and 
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dynamical evolution of galaxies in general. Understanding the birth 
and life of DGs on the other hand could prove crucial when trying 
to understand the formation and evolution of more massive stel- 
lar systems. And even more: cosmological simulations of galaxy 
formation exhibit very strong dependence on the way star forma- 
tion, feedback (FB) through supernovae a nd stellar winds (SW) 
are modelled within ste llar systems (see e.g. lWhite & Frenklll99ll ; 
iKauffmann et al.lll999l) . Therefore realistic simulations of single 
dwarf galaxies can supply information on how to implement star 
formation and feedback proce sses in large cosmological simula- 
tions (e.g. lSpringel et al.l2005T) . where a fully consistent implemen- 
tation of SF and FB should be possible for moderate simulation 
volumes. 

From the theoretical point of view various approaches for the 
study of DGs are possible, most notably there are semi-analytical 
models (SAMs), consisting of a limited set of analytical equations 
solved numerically, and fully numerical simulations, Lagrangian 
(e.g. N-body/SPH) or Eulerian (grid codes). SAMs have the ad- 
vantage of computational efficiency, due to their speed it is possible 
to scan a large part of the parameter space in a reasonable period 
of time. In numerical simulations however a fully self consistent 
implementation of SF, FB, chemical enrichment, gas cooling, ... is 
possible, at the price of high computational effort. 

iDekel & Silkl dl986h were the first to model Dwarf Galaxies 
(DGs), emphasising t he impact of supe r nova-d riven winds on the 
formation of dwarfs. lYoshii & Arimotd ([1987) found good agree- 
ment between observational properties of DGs and those predicted 
by their SAM. These authors calculated the properties of DGs (and 
globular clusters), using a model that takes into account the dynam- 
ical response to a supernova-driven wind. An evolutionary popula- 
tion synthesis method is employed to bring the models into the ob- 
servational domain. According to these models, spheroidal galaxies 
are a one-parameter family of birth m ass (initial total mas s M t .j) . 
These models were later extended in Nagas hima & Yoshiil 12004T) . 
who used a Monte Carlo technique t o construct the merger tree 
leading up to the formation of a DG. iNagashima et al.l 120050 in- 
corporated N-body simulations instead of a Monte Carlo method 
to construct the merger tree. They found good agreement between 
galaxies produced by their SAM, ranging from dwarf elliptical to 
elliptical galaxies, and observations. SAMs are also used to investi- 
gate certain processes involved in the birth and life of dwarf galax- 
ie s rather than trying to reproduce dwarf galaxies, e.g. the model 
of iBurkert & Ruiz-Lapuentd 11997|). investigating the role of SNIa 
in delaying star formation and of Ferra ra & Tolstovl (120000 who in- 
vestigated the role of stellar feedback and dark matter. 

Although plenty of numerical research has already been car- 
ried out concerning the formation and evolution of massive el- 
liptical galaxies (EG) (to name a few: |Ka tz 1991; Kawayj : 



Chiosi & Carraro 2002; Kawata & Gibson 2003; Kobavashi 2005; 



2001; 



Merlin & C hiosi 2006), starting from various cosmological initial 



conditions, only few people have tried directly modelling DGs 
with fully numerical simulations. One problem is that DGs are 
too small to be adequately resolved in large cosmological sim- 
ula tions. A first fully nume rical study of dwarfs was performed 
by iMori et al.l 1 19971 1 19990 . w ho explored the role of gas cool- 
ing in dwarf galaxy formation. IChiosi & Carrard 12002J) used an 
N-body/SPH code to simulate the formation of galaxies, includ- 
ing dwarf spheroidal, dwarf elliptical and normal elliptical galax- 
ies. They started from cosmologically inspired initial conditions. 
Assuming that star formation starts after an initial merging period, 
the simulations start with the colla pse of an initia lly homogeneous 
gas cloud in a dark matter halo IWhite & Reesl 1 19780 . Observ- 



able propertie s of DGs were broa dly reproduced by their models. 
More recentlv. lStinson et alj J2007T) conducted simulations of dwarf 
spheroi dal galaxi es (dSph). They used initial conditions similar to 
those of IChiosi & Carr aro ( 2002) , except that a deg ree of rotation 
was added to the gas. Like lChiosi & Carrar o (2002), thes e auth ors 
find burst-like star for mation. The work of IStinson et al.l 120070 is 
based on the work by IStinson et al.l 120060 . who investigated the 
influence of star formation criteria on the formation of a galaxy, as 
well as th e dependence of th e proce sses involved on the number of 
particles. iMashchenko et alj (2007J) performed fully cosmological 
simulations of the formation of a few dwarf galaxies, focusing on 
the gravitational heating of the dark matter, induced by gas flows. 
They argue that this heating can effectively convert a cusped dark 
matter halo into a D M halo with a central co re. Other studies of 
dSphs include those of iRead & Gilmord 120050 . who looked at the 
effect of baryonic mass loss from a two c omponent galaxy (dark 
and baryonic matter) and lRead et alj ([2006), who used cosmologi- 
cal simulations to investigate the smalles t baryo nic building blocks 
for galaxy formation. Mashchenko et alj 120050 employed a simple 
model for the formation of dwarf spheroidal galaxies. A single star- 
burst was imposed on a gas distribution, employing a density cutoff, 
after which the remaining gas was removed. They find reasonable 
agreement with data from local group dwarf spheroidal galaxies. 

A stu dy of dwarf galaxies in a cosmological context was per- 
formed bv lScannapieco et alj J200ll) , investigating the effect of gas 
outflows on the formation of dwarf galaxies. Other studies us- 
ing simulations, not directly aimed at r e produ cing dwarf galax- 
ies, include those of iMac Low & Ferrara d 19990 . who used a two- 
dimensional grid code to mod el blowout/blow - away o f gas in dwarf 
galaxies, and more recently iMarcolini et alj 120060 . who used a 
three-dimensional grid code to study the ISM in dwarf spheroidal 
galaxies by imposing a burst-like Star formation history (SFH) and 
comparing the results to the Draco dSph. 

In this paper we present the results of our models of dwarf 
galaxies, with initial conditions taken from the ACDM formalism. 
As we construct spheroidal, non-rotating galaxies, we expect them 
to resemble dwarf spheroidal and dwarf elliptical galaxies. We took 
great care to avoid the introduction of a direct dependence on the 
number of particles or on the resolution into the formalism. We also 
investigate the role of the different star formation criteria in detail. 
We then compare our model data to observations using the central 
velocity dispersion, the half-light radius, B— V broad band colours, 
total metallicity and the fundamental plane. 

The outline of this paper is as follows. In §[2] we summarise 
the most important features of the modelling code. The initial con- 
ditions of our models are discussed in § [3] In § [4] we compare our 
models with observational data. We end with a summary and our 
conclusions in §|5] 



2 THE CODE 

We cons tructed our models using the N-body/ SPH code 
HYDRA dCouchmanet~aIlll995l ; iThacker et al1l2000Q . In an N- 
body/SPH code the gravitational equations are solved with an 
N-body integrator. Hydrodynamical forces are applie d through a 
Smoothed Particle Hydrody namics (SPH) formalism ILucvll 1977c 
iGingold & Monaghanlll977l) . In this formalism the density at the 
position of the i-th (gas) particle is: 



Pe,i =^2mjW(nj,hi 



(1) 
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i | , hi is the smoothing length of the i-th particle 
and W is the smoothing kernel. The use of the SPH formalism lies 
in the fact that derivatives of e.g. p g can be expressed as derivatives 
of the smoothing kernel W. HYDRA was modified by us to include 
star formation, feedback, chemical enrichment, radiative cooling 
and supemovae (type II and la). 



2.1 Star formation 

To implement star formation in our code we use a phenomenologi- 
cal, Lagrangian approach. The first step is to formulate a set of star 
formation criteria. 



the gas density is higher, this check is performed more in highly 
clustered regions, which makes perfect sense when considering our 
SFC (eq.[2}{4}. This does however not imply a dependence of the 
particle number, because the newly formed star particles have a 
mass proportional to their parent gas particles. When e.g. doubling 
the number of gas particles (in a certain area, or in the entire sim- 
ulation), the mass attributed to each gas particle is cut in half. As 
our SFC are independent of the number of particles, two times as 
much gas particles will form stars, leaving the total mass of newly 
formed stars invariant. The number of gas particles in a certain area 
thus determines the average mass of star particles spawned in that 
area and the resolution with which the SFC are checked, it does not 
(or very lightly) influence the total star mass generated in that area. 



2.1.1 Star Formation Criteria 

We strive to select a set of Star-Formation Criteria (SFC) that al- 
lows our simulations to mimic re al-life star form ation. A number 
of implementations is compared in lKav et al . (2002), who find that 
different implemen tations of SF g ive comparable results. We will 
use a variant of the lKatzet alj 119961) criteria, where it is assumed 
that star formation takes place in cool, dense, converging and grav- 
itationally unsta ble molecular clouds. A good review of these cri- 
teria is given by IStinson et alj J2006I) . We select cool, dense and 
converging clouds in a very straightforward manner: the tempera- 
ture of gas eligible for star formation must be lower than a certai n 
critical temperature T c = 15 000 K l lStinson et al]|200d I2007T) . 
The density of a star forming gas cloud must exceed a critical den- 
sity p s = 2 x 10~ 25 g cm-3, n H « 0.1 cm - 3 JKawatall200ll ; 
IStinson et alj|2007l) . Furthermore the local gas flow must be con- 
verging: V • v < 0. 

An important feature of a set of SFC is that they should be 
independent of the number of particles. Not only is it a necessary 
condition that simulations converge, we also want them to converge 
as fast as possible to avoid extreme computation times. It is fur- 
thermore impossible to keep the mass resolution fixed in simula- 
tions spanning a wide mass-range (e.g. factor 100), because that 
would imply a proportional variation of the number of particles 
(e.g. 3 x 10 4 — 3 x 10 6 ), resulting in drastic changes in and large 
requirements for the number of particles. This reasoning rules out 
the commonly used implementation of the Jeans-criterion for dy- 
namical instability : td < t s , with td the dynamical time and t s 
the sound crossing time, defined as h/c s — t 3 . Here h is the SPH 
smoothing length, c s is the sound speed. Using h introduces a di- 
rect relation to the number of particles. One should moreover not 
attribute any physical meaning to it because it is an artificial con- 
struct. Furthermore, as td ~ l/\/pg and c B ~ yT, with p g the 
gas density and T the gas temperature, the Jeans-criterion is al- 
ready contained in our SFC. Therefore, instead of trying to find a 
meaningful length scale for the gas, we did not include a separate 
criterion for gravitational instability of the gas. Our SFC are thus: 



V- v < 

Pg > pc = 0.1 cm~ 3 
T To = 15000 K. 



(2) 
(3) 

(4) 



Each time step the SFC are checked for each gas particle. When a 
gas particle meets all three conditions it becomes eligible for star 
formation. Although one should refrain from attributing particle 
properties to the SPH gas particles because they are actually (only) 
co-moving grid points, this course of action is justified when inter- 
preted as a particle number independent check of the SFC through- 
out our gas cloud. As there are more grid points (particles) where 



2.1.2 Star Formation Law 

When we have selected those particles eligible for star formation, 
we have to prescribe a recipe to turn them into stars. For our SF law 
we adopt the Schmidt law of (1959) : 

dps _ _dpg Pg 

di di *i g ' [ ' 

where p s and p g are the density of stars and gas respectively, c* is 
a dimensionless scaling constant. i g Is a characteristic time-scale 
for the gas, we adopt i g = id, where id is the dynamical time. We 
define the dynamical time as: 



t d = 



V 47vG Pz 



(6) 



where G is the gravitational constant. Throughout our simulations 
we set c* = 1. It is possible to constrain c* and p c by fittin g 
star formation rates to the Kennicut-Schmidt law ( Ken nicuta 1998T) . 
IStinson et alj J2006I) found little influence of c* on the median SFR 
when going from 0.05 to 1 (their fig. 14). They opted to use 0.05, 
we use the canonical value 1 . 

We argue that th e use of p to t instead of p g , as argued in 
iBuonomo et alj |2000h. in the definit ion of the dynamical time is 
less justified. IBuonomo et alj ( 12000]) argue that the use of ptot 
results from deriving the Jeans criterion for instability (see e.g. 
iBinnev & Tremai ne 1987) for a gas cloud embedded in a dark halo. 
This is only true if the dark matter collapses along with the gas but 
not if the dark matter merely acts as a background fluid with, lo- 
cally, constant density. We therefore use only the gas density p g in 
our definition of id. 

Integrating equation (O over the interval At (assuming con- 
stant td), multiplying with a volume V and dividing by the initial 
gas mass M gj o gives the probability that a gas particle forms a star 
particle: 



exp 



At 

17 



(7) 



This equation has the immediate advantage that it introduces a 
quasi-independence of the number of time-steps into the formal- 
ism. If we bridge a time-interval Ai with one step, the P to form 
a star particle is: P = 1 — exp(— At/td), For small values of 
At/id we have: P « At/td- For the total probability to form a star 
particle using 2 steps we then have: (1 — Ai/(2id))Ai/(2id) + 
Ai/(2i d )(l - Ai/(2i d )) = Ai/i d - 2(Ai/(2i d )) 2 . This is, to 
the first order in Ai/id, equal to the single-step value. As Ai/id 
is typically 0.01 (see further), the approximation is justified. Equa- 
tion {7} is implemented using a Monte-Carlo procedure: a random 
number r between and 1 is drawn, if r < p* the gas particle 



4 S. Valcke et al. 



gives birth to a star particle. The mass of a newly formed star par- 
ticle is fixed at 1/3 of the mass of the parent gas particle. Each gas 
particle is allowed to form 4 star particles after which its remain- 
ing mass, metals, energy, ... ar e distributed amo ng neighbour gas 
particles, using the SPH kernel. iKennicuttl d 19981) finds that the me- 
dian rate of gas consumption is 30 per cent per 10 8 yr. Calculating 
the average central dynamical time td,c (equation ® for gas par- 
ticles within a radius of 0.5 kpc from the centre) we find a value 
O(10 8 yr). The formula for the average mass of a gas particle after 
N equal time-steps, taking into account the efficiency of 1/3 for SF 
as well as equation ifTJ, is: 



M = M 



i+£ 



min (i, 4) 



- 1 



p*(i-p* 



0]\ 



,(8) 



where & N = N\/(i\(N - i!)), p*(l - p*) N -*C l N is the probabil- 
ity that a star particle will form stars i times in N steps and (2/3)' 
stems from the 1/3 efficiency of our SF. When using equation {8) 
we assume that the gas particle satisfies the 3 SFC for the dura- 
tion of the N time-steps. We also assume that p* is invariant during 
these time-steps. The limited number of star formation episodes (4) 
is taken into account by using min(i, 4) instead of i, which comes 
down to setting the efficiency in terms of mass of more than 4 SF 
episodes equal to 4. This yields a negligible correction. We can 
write equation as p* ~ 1 — exp(— 1/N), with N = id/At. As 
td 0(1O 8 ), and the time-steps in our simulations are 0(1O 6 ), we 
have N 0(100). Using this value gives M « 0.72 M , so our fac- 
tor 1/3 models 30 per cent efficiency rather well. This is of course a 
crude estimation, but we preferred it over adding a parameter to our 
model, which would multiply the required number of simulations 
by a further factor. 

2.2 Feedback 

There are two mechanisms through which stars can return met- 
als to their environment: stellar winds and supernovae. These can 
be modelled given knowledge of the initial mass function (IMF), 
which determines the distribution over mass within a "star particle" 
(single age single metallicity particle or SSP). We have adopted the 
Salpeter IMF in our simulations: 

*(m)dm = Am' (1+X) dm, (9) 

where x — 1.35 and A is fixed by the chosen normalisation. 
<E>(m)(= diV/dm) is the probability that a star with mass m 
resides in the SSP. We normalise the total probability to one: 



/™ u $(m)dm 
we have A — 



1. Using mi = 0.1 Mq and m u 



60 M 



© 



0.06. To calculate the average value of a certain 
quantity B of a SSP we then have: 



<£> = 



B(m)$(m)dm x 



Mssf 



/*"" m$(m)dm' 



(10) 



with ttib.i and mp,, u respectively the lower and upper bound for the 
mass interval where B applies, and Mgsp the mass of the stellar 
particle. Application of equation dlOt for energy feedback of SN II, 
where we assume B(m) — -Bsnii to be independent of mass, and 
with m S Nii,i = 8 Mq, m S Nii,u = 60 Mq, gives: 



E u 



Est- 



l e 



$(m)dm 



a © 



M SS p 



' 60 Mq ' 

io.i mq m$(m)dm 



(11) 



which gives us: 

Stot.smi = 7.31 x IO^MsspSsniiMq 1 



© 



(12) 



For SWs the derivati on is, apart from repla cing -Esnii by i?sw, 
identical. Following iThornton et alj dl998T) we set Ssnii = 
10 51 erg, -Esw is set to = 10 50 erg. The actual energy re- 
turned to the ISM, for SNe as well as for SW, is implemented as 
£f b x Etot, where £fb was chosen to be 0.1, following the results 
of IThornton et alj dl998l) . The returned mass fraction Jsnii can be 
easily calculated by subtracting the total mass that is not returned 
to the interstellar medium (ISM) by exploding stars from the to- 
tal mass of stars eligible to go supernova. With the approximation 
that the mass remaining in dark objects (black holes, neutron stars, 
white dwarfs) after a star goes SN II is constant, M rom « 1.4 Mq, 
we have: 



-Fsnii = 



pSNII.u m $( m)dm . 
Jm SNII,l V ' 



M r< 



, f msNII ' u *(ra)dm 

Jm SNII,l V ' 



f u m$(m)dm 



= 0.112. 



(13) 



The yield of element i by SN II is calculated as (equation ( 1 101 

r sNII '°M,(m)$(m)dm 

J ™SNII,1 V ' y ' 



Mi = Mssp- 



J m " m&(m)dm 



(14) 



For the metal yields of stars w ith mass m we fitted c ubic splines 
to the data points tabulated bv lTsuiimoto et al.1 (1995J). For SN la 
the structure of the used formulae is analogous, the only difference 
being the inclusion of a factor AsNia in every right hand side of 
the above equations. AsNia represents the fraction of stars in the 
mass range of SN la that actually go supernova, because SN la are 
believed to take place after a period of Roche lobe overflow in a 
binary star system with a white dwarf and a red giant companion. 
We calculate ^sNia as follows (msi-jia, the lower limit for SN la, 
is taken to be 3 Mq, the upper limit for SN la is msNii.i): 



A^SNIa 



= 0.15 = As* 



r-sNii.i $( m ) dm 



(15) 



Ar SN „ """"*f msNII ' u $(m)dm' 

J ™-SNII,1 V ' 

where 0.15 was found by iTsuiimoto et al.1 d 19951) as the number 
ratio that best reproduces the observed abundance pattern among 
heavy elements in the solar neighb ourhoo d. This gives us: ylsNia = 
0.0508 (compare to e.g. iKawatd J200lh who use 0.04). We note 
that because SN la take place when white dwarfs reach the Chan- 
drasekhar limit (1.4 Mq), equation dl4t will be simplified for SN 
la because Mj(to) becomes Mi, independent of mass. 

Using AsNia and an equation analogous to equation d 1 3fc . set- 
ting Mrem = 0, we find: FsNla = 0.00502. Using the analogon of 
equation d 14b for SN la, keeping in mind that the yields per super- 
nova are independent of the progenitor mass, we find: Mj,sNla = 
Fj,sNiaMssp with Fj,sNia = O.OO llmi.sNiaMg 1 . We used the 
yields given in lTravaglio et alj (12004) (their b20_3d_768 model). 

For the time intervals for SNe and stellar wind we use a 
formula for the m ain-sequence lifetime of a star with mass m 
JDavidetalJl990l) : 



log r(m) = 10 - 3.42 log m + 0.88 (log(m)) 2 



(16) 



Filling in the various lower- and upper bounds for the mass of SNe 
we have: SW: 0-4.3 x 10 7 yr, SNe II: 5.4 x 10 6 -4.3 x 10 7 yr, SNe 
la: 1.543 x 10 9 - 1.87 x 10 9 yr. For SN la a net dela y of 1.5 Gyr was 
applied, following the result of lYoshii et alj ( 1199a) . who found 1.5 
Gyr for the mean lifetime of SNe la based on a chemical evolution 
model applied to the solar neighbourhood. This time interval is the 
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result of the time needed for the accretion of mass by the WD to 
push its mass up to the Chandrasekhar limit. For the distribution 
of feedback energy we used the simplest prescription: a constant 
amount of feedback over the time intervals. 

The actual feedback uses the SPH smoothing kernel: if star 
particle i with smoothing length h s distributes an amount of energy 
E, then a neighbouring gas particle j receives an amount of energy 
AEr. 



AE, 



rrijWQf'i -fj\,h s )E 



E: =1 rn h W(\ 



rt-rk 



,h a 



(17) 



Metallicity-dependent cooling (with a maxi mum metallicity of 1 
Z^)i s implemented using the tables given bv lSutherland & Dopital 
( 1993), allowing our gas to cool to a minimum temperature of 10 4 
K (cooling beyond this temperature is possible however through 
adiabatic expansion). To compensate for the fact that we are un- 
able to resolve the hot, low density cavities in the gas, origi- 
nating from SN explosions, we implemented a period where ra- 
diative cooling is switched off, only allowing adiabatic cooling 
JThacker & Couchmanll200d) . To be precise a particle is not al- 
lowed to cool radiatively during a certain time-step if it was heated 
by SNe during that ti me- step. 

As explained in iBate & Burkern 019971) . a difference of the 
gravitational softening (e) and SPH smoothing (h) can lead to arti- 
ficial convergence/divergence of SPH particles. We tested e — 30 
pc and e = 60 pc and found very comparable simulation results, 
permitting us to use 60 pc. h ranges from on average 120 pc (C02 
model) to 60 pc (C09 model), so e and h are of comparable size. As 
the global evolution of our models is determined by the in-fall into 
the deep potential well, we do not have to be very concerned about 
mismatching of e and h on that scale. And as the resolution of our 
simulations is insufficient to capture the collapse of small (O(30 
pc)) molecular clouds, we need not be concerned about artificial 
convergence/divergence of SPH particles on that scale. 

Broad band colours of our models a re calculated (with bi - 
linear interpolation) using the models of IVazdekis et alj (1996), 
who provide Mass/Luminosity values for SSPs according to metal- 
licity and age. 



3 INITIAL CONDITIONS 

As initial conditions, we use a relatively simple setup, assuming a 
flat A-dominated cold dark matter cosmological model. We fur- 
ther assume that in the process of dwarf galaxy formation, star 
formation comes into play only after there has already been sub- 
stantial merging. Support for this assumption can be gleaned from 
Nagashima & Yoshii (2004), who find that SAMs with long star- 
formation time-scales at high redshift provide the best agreement 
with the observed properties of dwarf galaxies. Furthermore, we 
used th e freely available res ults from the milli-Millennium Simu- 
latiorj 1 ] ( Springel et al. 2005) to assess the formation redshift of the 
halos that at redshift z = fit the description of an isolated dwarf 
galaxy halo. The dwarf galaxies modeled here are actually below 
the resolution of the Millennium Simulation, which requires a min- 
imum of 7V P = 20 particles, or a mass of about 1.7 x 10 10 Mq/i _1 , 
to properly identify a halo. This corresponds to the mass of the C09 
halo, our most massive model (see table[3}). 

We queried the milli-Millennium Simulation for the formation 
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Figure 1. The number of halos drawn from the Millenium Simulation with 
particle numbers a) 175 < N p < 185 and b) 35 < N p < 45 at z = 0, 
binned as a function of how many gigayears ago these halos first reached 
particles number above a) 100 and b) 20. Halos that undergo a major merger 
after their formation are explicitly excluded from this tally. 



redshift of halos that at z — would be associated with isolated 
dwarf galaxies. We defined the formation redshift as the redshift 
at which a halo achieved approximately 2/3 of its final mass, al- 
lowing for a factor of 1.75 of further mass growth through minor 
mergers. Halos that undergo a major merger after their formation 
are explicitly excluded from this tally. We cannot derive the forma- 
tion redshifts of halos with particle numbers N p ~ 20 at z — 
because (i) the closer the halo particle number gets to the lower 
limit of halo detection, the less accurate we can determine the ac- 
tual time of birth of that halo and (ii) for such low particle num- 
bers it is impossible to construct a proper merger tree. Instead, we 
estimate the formation redshifts of halos with particle numbers be- 
tween 35 and 45 at z — by querying when they first encompassed 
at least 20 particles. In order to investigate whether the halo forma- 
tion redshift depends on galaxy mass, we also queried the milli- 
Millennium Simulation for the redshift at which halos with particle 
numbers 175 < N p < 185 at z — first encompassed 100 par- 
ticles. This again alows for a factor 1.75 of further mass growth 
through minor mergers. 

It is clear from Fig.Q]that the peak of halo formation is situ- 
ated at about 10 Gyr ago, with the majority of halos having formed 
between 8 and 1 1 Gyr ago. A closer inspection of both panels of 
Fig. Q] reveals that the formation rate of the least massive halos, 
presented in panel b), peaks at ~ 10 Gyr ago. This is about 1 Gyr 
earlier than the formation rate of their more massive analogs, pre- 
sented in panel a). Extrapolating towards the even lower masses of 
the galaxies we modeled for the present paper, one could expect 
the halo formation time to shift to even earlier times, with the old- 
est halos being formed as early as 12 Gyr ago. We will, therefore, 
assume that the dark-matter halos of the dwarf galaxies we mod- 
eled formed at a redshift z = 4.3 and neglect further growth of the 
halo. Furthermore we neglect any environmental effects (e.g. ram- 
pressure stripping) on the chemical and dynamical evolution of the 
models. These environmental in fluences are not always negligible 
(see e.g. lPenarrubia et alj ( 120070 ). Thus, the models describe dwarf 
galaxy formation as the isolated collapse of a gas cloud in the po- 
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tential well g enerated by a spheric al dark matter halo and the gas 



tential well generated by a sphenca 
itself (see e.g. lGao & Theunsll 20071) 



We use the following cosmological parameters: h = 0.71, 

foot = l.pmatt er = ^m = 0.2383, S7dM = 0.1967 

dSpergel et alj|2007l) . where h is the normalised Hubble constant. 
The mean density of the Universe as a function of redshift z and 
Hubble constant Ho = 100 h is: 

3/i 2 100 2 r ^ ,„ 
P u = _^ [S2 m (l 4 



SttG 



= 1.99 x 10~ 2 V [fi m (l + z) 3 



- I - On] 

l-o, 



gem 



(18) 



We take the over-density of matter to p u to be 5.55, the value when 
the local flow detaches itself from the Hubble flow, according to 
the Tolman model of a spherical over-density, z at the start of the 
simulations is taken to be 4.3, corresponding to approximately 1.5 
Gyr after the birth of the universe. The simulation length is set to 
10 Gyr, which is long enough to cover the most important part of 
the star forming period of the simulated galaxies. The density of 
matter at the start of the simulations is then: 



pm = 5.55 x p u . 



(19) 



When the total mass and its density are known we can fix the outer 
radius: 



Rs.&&(z) — 



1/3 



M T 

Pm 



1/3 



0.09617 
5.55V3 



Mi 



-,1/3 



h 2 (On(l + Zf + 1 - On) 



(20) 



where Mr is expressed in Mq and -R5.55 in kpc. At the start of the 
simulations the gas particles are randomly distributed in a sphere 
with radius given by equation ( 1201 ). The gas particles are initially 
at rest, their initial metallicity is set to 10 -4 Zq and their initial 
temperature is 10 4 K. 

Several choices are possible for a DM pro file. There is 
the N FW profile, deduced from simulations in iNavarro et al.l 
dl997l) . However observations as well as simulations seem to 
imply that the DM halo has a central core rather than a cusp 
dSalucci & Persidl997l ; lMerritt et alj20o3 : lGentile et alj200l and 
references therein). Oth er previously used profiles include the mod- 
ified isothermal sphere (Binney & Tremaine 1987) or the King pro- 
file (e.g. iMori et all 17999), see Binney & Tremaine ( 1983). For 
the DM halo we choose the Kuz'min Ku tuzov (KK) profile as 
presented by iDeionghe & de Zeeuwl (j_1988). The KK profile has 
several distinct advantages : it is, apart from a numerically well- 
behaved integral, analytical, and it has a central core. Furthermore, 
Sellwood & Valluri ( 1997) have carried out a stability analysis of 
the KK halo, concluding that haloes with flattening up to E7 are 
stable. The density of this model is given by: 



Pkk{R,z) = 



Mc 2 (a 2 + c 2 )R 2 + 2a V + 
4tt (a 2 c 2 + c 2 R 2 + a 2 z 2 ) 3 / 2 



2a 2 c 2 + a 4 + 3a 2 Va 2 c 2 +c 2 R 2 +a 2 z 2 



(R 2 + z 2 +a 2 + 2Va 2 c 2 +c 2 R 2 +a 2 z 2 ) 3 / 2 



(21) 



with a and c parameters controlling the length of the semi-major 
and semi-minor axis respectively, R is the distance in the a;j/-plane 
and z is the distance normal to this plane. 

We generated the DM halo using a standard Monte Carlo 
sampling technique. For each particle we first generate three co- 
ordinates R, $, z, distributed as pKK (equation d21t ). Then vr, v$ 
and v z are found using a standard acceptance-rejectance technique 
using the distribution function of the KK model. The halo den- 



sity is set to zero outside the equipotential surface $kk(-R, z) = 
*KK(^cut,0), with Rcut = 25 (units (a + c)). This large value 
ensures that the effect of this cut-off on the halo's stability is negli- 
gible. We explicitly checked the halo's stability by letting it evolve 
in isolation for the same duration as the science simulations. 

Now that we have fixed the cosmological parameters (includ- 
ing the ratio of DM mass to gas mass), we are left with 3 parameters 
determining the model: Mtot , a and c. c determines the amount of 
flattening we want for the halo, where a — c gives a spherical halo. 
As we restrict ourselves here to spherical systems we set c = a. To 
eliminate a further parameter we set, analogous to -R5.55 in equa- 
tion J201. 



cVm, 



kpc 



(22) 



where C is fixed by associating one value of a with a value of Mtot . 
Several values of C are explored with our models. This leaves us 
with one parameter to vary in a series of simulations: M to t- 



4 RESULTS 

In the following sections we discuss the models. Information about 
these models can be found in Tables Q] (A models) and [3] (C and 
D models). The A models all have an identical setup, apart from 
the number of gas particles. The number of DM particles is set to 
15 000. We use these A models to determine the minimum number 
of gas particles that is required for the simulations to be numerically 
converged. The C models represent the "real" models. We com- 
pare in detail some structural parameters of the simulated C-model 
galaxies with those of observed dwarf galaxies, in order to assess 
how well the simulations reproduce reality. The B models are es- 
sentially a larger version of the C models, being identical in setup, 
apart from a smaller value of C (equation (22): C = 0.000928. 
We also sampled the mass range differently: the total mass of the 
Bx model does not necessarily correspond to the mass of the Cx 
model (B models: 0.25, 0.3, 0.65, 1, 3, 6 xlO 9 Mq). This how- 
ever has no influence on a comparison of the B and C models. The 
B models are included only to illustrate the effect of varying C on 
the half-light radius and the central velocity dispersion of the model 
galaxies (sections l4.4l and l4~6l . They were constructed using 15 000 
dark matter particles. The two D models do not obey equation ( 122b . 
They were constructed using an arbitrary (high) value of C. The 
D models are included to show that it is difficult to extend our C 
models beyond the mass where they break down by fiddling with 
some parameters. 



4.1 Number of particles 

To be able to compare simulations covering a wide range of initial 
masses, the simulations have to converge numerically. At a certain 
point a further increase of the number of particles should make no 
(or very little) difference. For this reason we took special care in 
selecting the SFC (equations [2}|4}. Results of a first set of simula- 
tions for a range of particle numbers (A models) are shown given 
in Figs.[2}{3] These models all use the parameters of the C05 model 
(see table |3). We use the models listed in table Q] to investigate 
the dependence of the models on the gas-particle number. These 
models are all different realisations of the same set-up: that of the 
C05 model (see table [3)- The 30k(rd), 50k and 75k sets consist of 
5 simulations; the 30k set of 6 simulations and the 100k set of 3 
simulations. The simulations of the 30k set all start from the same 
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Table 1. Information concerning the simulations with varying particle num- 
bers. See text for more details. 



Name 



iY„ 



Comments 



30k 


3.0 x 10 4 


30k(rd) 


3.0 x 10 4 


50k 


5.0 x 10 4 


75k 


7.5 x 10 4 


100k 


10.0 x 10 4 



Identical initial conditions files (ICfs). 

Resampled ICfs (different Poisson noise). 

Identical initial conditions files. 

Identical initial conditions files. 

Identical initial conditions files. 



Table 2. Information concerning the simulations with varying dark matter 
particle numbers. (-/V gas = 30 000) 



Name 



^DM 



# runs 



15kdm 


1.5 x 10 4 


5 


30kdm 


3.0 x 10 4 


5 


5k-to-25k 


(0.5,1,1.5,2,2.5) x 10 4 




30k-to-60k 


(3,3.5,4,4.5,5,6) x 10 4 





0.006 




Figure 2. Star Formation Rate (SFR, in Mq yr _1 ) as a function of time 
(Gyr) for simulations with varying particle numbers. See text and table [T] 
for more details. 



initial conditions file, only the random number seed of the stochas- 
tic star formation recipy (eq. 0) varies between simulations. The 
same goes for the 50k, 75k and 100k sets of simulations. For the 
30k(rd) set, the initial conditions file was resampled for each simu- 
lation, i.e. they have variable Poisson noise at constant gas-particle 
number. The number of dark matter particles was kept constant at 
30 000 particles for all simulations. 

A first indication of the global behaviour of the models is 
shown in Fig. [2] where we plotted the mean SFRs of the sets of 
models listed inQ] with la "errorbars" quantifying the variation 
between the simulations belonging to each set. As a reference, the 
SFR of the C05 model is also plotted. Obviously, the global SFRs 
are qualitatively and quantitatively in good agreement and they do 
not change significantly with gas-particle number. Of course, slight 
differences between the different sets of models during the first 
star-formation peaks are amplified towards the later peaks. In fig- 
ure [3] we plot different observables as a function of stellar mass 
for the models listed in table Q] as well as the C05 model. We also 
plot a line connecting the C-series models. Apparently, simulations 
with higher numbers of gas particles tend to form slightly more 
stars, the difference between the 30k and 75k sets of simulations 
being of the order of 10% in terms of stellar mass and with com- 
parable differences on the other observables. There is however no 
straightforward converging behaviour as the 100k set forms on av- 
erage less stars than the 75k set. The least well constrained quan- 
tity is the percentage of retained gas, with a 30% increase between 
the 30k and 75k models. Clearly, the scatter between the individ- 
ual models diminishes towards larger gas-particle numbers but in- 
creasing the gas particle number beyond 75 000 seems to have no 
further effect. We conclude that if one wants to compare an indi- 



vidual model in detail to an observed galaxy, it seems advisable to 
use at least 100 000 gas particles in order to minimize the model- 
to-model scactter. For a comparison between observed and model 
photometric and kinematic scaling relations, a gas-particle number 
of 30 000 should suffice. The model-to-model scatter is sufficiently 
small and, in any case, this scatter moves models along the scaling 
relations so the latter will not be significant ly affected. This conclu- 
sion is a lot stricter than the one reached in lLia et alj I l2000h . where 
the requirement of comparable star formation histories resulted in 
iV gas > 10 000. 

To assess the number of dark matter particles to be used in the 
simulations, we constructed a further series of models (see table[2]l. 
For the 15kdm we started from the same initial conditions file and 
used 5 different random number seeds for the star formation recipy. 
The same was done for the 30kdm models. For the 5k-to-25k and 
30k-to-60k models the number of dark matter particles was varied 
whilst retaining the initial conditions for the gas particles. Only 1 
run per initial conditions file was done for these models. Results 
are shown in Fig. [4] Two important conclusions can be drawn from 
this figure: (i) there is no significant decrease of the scatter on the 
individual models as a result of an increase in dark matter particles 
and (ii) DM particle numbers of 15 000 and up seem to be sufficient 
to capture enough physics to closely follow the model relations. 
Indeed, apart from the leftmost two blue points in each plot (which 
correspond to 5 000 and 10 000 DM particles) all the models follow 
the C model relations. Although we do not rule out a decrease of the 
scatter on individual models as the number of dark matter particles 
is further increased, the above two results allow us to use 30 000 
dark matter particles to investigate the relations our model galaxies 
follow. 

Moreover, we performed a series of simulations with C09 ini- 
tial conditions but with the number of gas particles increased from 
30 000 to 40 000, 50 000, and 60 000. Due to star formation, the 
60k run reached a total of 200 000 particles after 5 Gyr and be- 
came impracticably slow and we decided to stop it. Comparing 
these models, we found that no observable scales directly with the 
number of gas particles and that the model-to-model scatter on the 
observables is too small to affect our results and conclusions. 

In conclusion: for our production run models (C and D mod- 
els) we use 30 000 gas particles and 30 000 dar k matter particles. 
These numbers used in are an improvement over lChiosi & Carrard 
(2002), who use 10 000 particles for dark and baryonic matter. 
They are less than t he 100 000 particles for DM and BM used by 
IStinson et alj ( I2007T) . but we have nevertheless demonstrated that 
this resolution is sufficient to capture the essentials of the physics. 



4.2 Surface brightness Profile 

The surface brightness (SB) profile of elliptical galaxies is de- 
scribed well by the de Vaucouleurs R 1 ' law, whereas dwarf el- 
liptical galaxies are fitted better by an exponential law (see e.g. 
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Figure 3. Variation of model observables as a function of gas particle number, shown for the models listed in table[T] The dash-dotted line connects the points 
of the C models. The symbol at the center of the error-bars indicates to which models these bars belong (it is not a datapoint). 
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Figure 4. Variation of model observables as a function of the dark matter particle number. The dash-dotted line connects the points of the C models. The 
symbol at the center of the error-bars indicates to which models these bars belong (it is not a datapoint). The leftmost two blue points in each figure are the 
models with 5 000 and 10 000 DM particles (from left to right). For clarity: the larger of the two sets of error-bars in each figure belong to the 30kg-30kdm 
models. See table|2]for more information about the models. 
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Table 3. Summary of various physical quantities for models C and D (at t = 10 Gyr). (1) Model name (2) Initial gas mass (IO 6 Mq) 
(3) Initial dark matter mass (IO 6 Mq) (4) a for the DM halo (kpc)(see eq. 12 U ) (5) Absolute magnitude in the blue band (6) Absolute 
magnitude in the visual band (7) B — R (8) Total final star mass (IO 6 Mq) (9) Percentage of initial gas mass remaining (< 30 kpc) 
(10) Percentage of gas blown out (> 30 kpc) (11) Half-light radius (kpc) (12) One dimensional central velocity dispersion (km s-1) (13) 
Metallicity (Z & = 0.02) 



Run M„ 



Mr 



Mu 



Mx 



B- R M B - 



%M, 



i BO R. 



°"1D,. 



Z(Z Q ) 



C01 


44 


206 


0.439 


-8.12 


-8.72 


1.067 


0.44 


21.0 


78.0 


0.16 


9.6 


0.048 


C02 


52 


248 


0.466 


-8.97 


-9.51 


0.981 


0.81 


19.5 


78.9 


0.18 


11.5 


0.052 


C03 


70 


330 


0.513 


-9.89 


-10.44 


0.997 


2.0 


15.6 


81.6 


0.23 


14.0 


0.045 


C04 


105 


495 


0.587 


-11.14 


-11.68 


0.973 


6.2 


9.1 


85.0 


0.31 


18.9 


0.073 


C05 


140 


660 


0.646 


-12.56 


-13.03 


0.878 


18 


17.9 


69.1 


0.48 


25.7 


0.174 


C06 


175 


825 


0.696 


-13.58 


-14.02 


0.831 


38 


32.7 


45.7 


0.67 


31.5 


0.230 


C07 


262 


1238 


0.797 


-14.20 


-14.83 


1.105 


122 


36.6 


16.5 


0.76 


39.3 


0.379 


COS 


349 


1651 


0.877 


-14.97 


-15.56 


1.043 


228 


22.7 


11.5 


0.62 


43.0 


0.535 


C09 


524 


2476 


1.004 


-15.11 


-15.77 


1.162 


393 


17.4 


6.80 


0.57 


46.7 


0.620 


D01 


873 


4127 


4.000 


-16.04 


-16.61 


1.023 


579 


29.0 


4.2 


1.13 


35.2 


0.470 


D02 


873 


4127 


6.000 


-16.37 


-16.82 


0.843 


488 


39.4 


4.4 


1.26 


30.9 


0.373 



Ijerien & Binggelilll9971 : lGraham & G uzman 2003, and references 
therein). Both profiles are encompassed by Sersics R 1 ' n law: 



I(R) = Io exp 



(l/n) 



(23) 



In Fig.[5]we present a plot of the surface brightness profiles of the 
C model galaxies as well as their respective best fitting Sersic pro- 
files. In Fig.[6]we compare the Sersic parameters n and /ip (or Jo ) 
for the models with d ata t aken from Graham & Guzman 120031) . 
|Pe Riicke et~ai] J2008h and iMieske et alj J20071) . As seen in Fig. 
[5] the correspondence with a Sersic profile is excellent. The more 
massive models (Mb < —13 mag) also have SB profiles that 
closely follow the Sersic law. Comparison with data in Fig.|6]for the 
C models shows that n as well as /io are in good correspondence 
with observations. A variant of the third Sersic parameter Ro, the 
half-light radius R c , is discussed in section l4l4l The D models per- 
form worse than the C models : whilst /io is only somewhat too low, 
the values for n are too low. This indicates that the SF should be 
more centrally concentrated. This is easily understood when con- 
sidering Table[3] by construction the D models have a large radius, 
naturally decreasing the resulting /io and n. 



4.3 Star Formation Rate 

The SFR of galaxies is an important tool for the study of their 
evolution. Observations indicate that the Star Formation His- 
tories of dwarf galaxies are varied and have a burst-like na- 
ture (e.g. ISmecker-Hane et al. 1996; Mated ll 9981 . and references 
therein). This phe nomenon was successful l y reproduced by N - 
body simulations JChiosi & Carrard 120021 : IStinson et al.1 2007). 
iDellenbusch et alj J20071) find ongoing star formation in their sam- 
ple of observed dwarf galaxies. They also note and an absence of 
any indication of a recent major interaction. These observations can 
be readily explained by invoking the burst-like nature of star forma- 
tion in dwarf galaxies. 

In Fig. [7] we plot the SFRs of the C models. We immedi- 
ately see that the SFR is of a burst-like nature. As there is 10 
to 30 per cent of the initial gas mass left at the end of the sim- 
ulations (Table [3} there is still fuel for further star formation. 
Most dwar f elliptical (dE) galaxies are found to be almost de- 
void of gas JYoung & Lcll997l : lMatecll998l : IConselice et al.l2003l : 




Figure 5. Surface brightness (SB, in mag arcsec ) as a function of radius 
(kpc) for simulations with varying initial mass. The full lines are the best 
fitting Sersic profiles. 



|Pe Riicke eta?] 120031 : iBouchard et"al] 120051 : iBuvle et alj 120051) . 
There is however a very stron g morphological segregation of dwarf 
elliptical galaxies dBinggeli et al.ll 19901) . These authors found that 
almost all dEs are found in clusters or in dense groups, or as com- 
panions (satellites) of a massive galaxy. This implies that dwarf 
elliptical galaxies, found in dense environments, are subject to sig- 
nificant ram-pressure s t ripping which deplete s their gas reservoirs 
JMori & Burkertll2000l : rMarcolini et alj|2003l) . This physical pro- 
cess is not included in our simulations so the modelled galaxies 
retain part of their gas reservoir. 

In Fig. [8] we plot the temperature distribution of the gas parti- 
cles as a function of their density, for the C05 model, at 4 different 
times: 1, 2.6 , 7 and 8.5 Gyr. In each figure the lower right quadrant 
contains the gas particles that satisfy the density and temperature 
criterion for star formation (eqs. and ©)• From Fig. [7] we can 
see that at 1 and 7 Gyr the galaxy is undergoing a star-burst. At 2.6 
and 8.5 Gyr there is no star-formation. From Fig.[8]it is clear that the 
quiescent intervals between star-bursts are caused by a decrease in 
gas density. This in turn is a result of the stellar feedback: the entire 
central population of gas particles is pushed to lower densities, in- 
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Figure 6. Upper panel: Sersic index n as a function of blue broad 
band magnitude M-q. Lower panel: central surface brightness fio (mag 
arcse c~ 2 ) from the S ersic fit to th e models as a fu nction of Mb . Data taken 
from Graham & Guzman 1 2003), De Rijcke et al. 12008) and (lower panel 
onlv HMieske et ah (2007) (the canonical value of 0.7 was added to absolute 
magnitudes and surface brightnesses in the V band to mimic their respective 
counterparts in the B band). 




Figure 7. SFR (Mq yr-1) as a function of time (Gyr) for the C models. 



hibiting further star formation. On the lower left part of each figure 
we can see gas particles being blown out of the galaxy, resulting in 
a declining density and temperature. On panels b), c) and d) we can 
clearly see the outward propagation of the gas particles expelled by 
the different starbursts towards the lower left corner of the figure. 

The cause of the burst-like behaviour can be seen on Fig. [9] 
where the effects of the various SFC are plotted for model C06. 
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Figure 8. The logarithm of the temperature, T, versus the density, p, for the 
gas particles, at 4 different times (C05 model). Particles are colour-coded 
according to their distance from the galaxy centre r. 



From the upper panel we learn that the temperature criterion (eq.O 
has little influence on SF, prohibiting a small fraction of the gas par- 
ticles from forming stars. The divergence criterion (eq. [5} is more 
restrictive, allowing on average 30 per cent of the gas particles to 
form stars and restricting the star formation to the more central re- 
gions of the galaxy. The real control over star formation resides 
however in the density criterion (eq. [4}. Not only is this criterion 
very restrictive (which of course depends on the value of p c ), typi- 
cally allowing 5-10 per cent of the gas particles to form stars, it is 
also responsible for the burst-like behaviour the star formation. Star 
formation is in a sense a self-regulating process: high density gives 
high SF, resulting in a lot of SNe which will heat up and disperse 
the gas, leading to the end of the SF episode. As the temperature 
criterion has only little effect it could be dropped altogether (no un- 
realistic star formation should take place in hot regions because hot 
clouds with high density will cool fast a nyway, hot clouds with low 
density are not allowed to form stars). iKravtsovl ( 120031) : lli et al.l 
J2005L 120061) use only the density criterion in their N-body/SPH 
simulations, whil e using enough particl es to satisfy the Jeans res- 
olution criterion (Bat e & Burkert 1997) and thus resolving Jeans 
instability in molecular clouds. They find that the Schmidt law (eq. 
[5} is a result of gravitational instability and thus of the overall den- 
sity distribution, which validates our SF because we also find it to 
be -mainly- determined by the density distribution. When compar- 
ing panels a) and b) of Fig. [9] we see the behaviour expected from 
td (equation l(6]l): peaks when there are few particles with high den- 
sity and vice verse, with an approximately constant value otherwise 
(« 0.12 Gyr). When going from panel 3) to panel 4), equation {7} 
is applied. Because of the moderate variations of td there is little 
qualitative difference between the two panels, although the quan- 
titative difference is large (factor 100). The use of equation Q is 
thus twofold: on the one hand it serves as a means to drastically 
reduce the number of star forming gas particles, on the other hand 
it introduces an independence of the number of time-steps into the 
star formation prescription. 

In Fig.[l0]we show a measure of the distribution of SF in space 
for our models. The apparent rise of the SFR towards larger radii 
(within a starburst) is caused by the rise in volume of the spheri- 
cal bins. The volume-averaged SFR, plotted in Fig.QT] shows that 
there is actually more SF in the centre of the Model galaxies. We 
can also see that in our more massive models, during the first star- 
burst, SF takes an increasing amount of time to reach the outskirts 
of the galaxy. On panels a), c) and d) of Fig.QJJ]we see that star for- 
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Figure 9. Data taken from the C06 model, a) Percentage of gas particles 
satisfying a SF criterion, b): Average dynamical time (Gyr) for gas particles 
within 0.5 kpc of the galaxy centre, c) Percentage of gas particles satisfying 
the 3 SFC. d) Percentage of gas particles forming stars. 



mation becomes more centrally concentrated as time passes: each 
subsequent star-burst is confined to a smaller radius. For model C02 
(panel d)) this effect is even stronger: the first star-burst contains 
two phases, with SF in the first phase extending to about fa 0.5 
kpc, in the second phase to w 0.3 kpc. Our mode ls thus naturally 
reproduce the observations of iTolstov et al.l ( 120040 . These authors 
found two distinct ancient components in the Sculptor dSPH, the 
metal-poor component being more extended in space. They con- 
sidered the interplay of gas infall due to gravitation and supernova 
feedback to be a possible explanation. iBattaglia et alj 112006) found 
a similar effect in the Fornax dSph: they found 3 distinct stellar 
populations (ages > 10 Gyr, 2-8 Gyr and < 100 Myr), where the 
younger, more metal-rich populations are found to occupy a smaller 
region of space. We can immediately qualitatively compare these 
findings with our C04 Model, which also gives rise to 3 different 
populations, with each subsequent population b eing more metal- 
rich a nd confined to a smaller region in space. iMashchenko et al.l 
(2007) offer a different explanation: the DM, gravitationally heated 
because of gas motion, in turn heats the stellar population. As the 
older (metal poor) population has somewhat more time to be heated 
than the metal rich population, it will occupy a larger region in 
space. These authors, however, do not provide quantitative predic- 
tions based on this scenario. Note that the apparent absence of SF 
in the inner radial bin, visible in FigQT| is an artifact caused by 
small number statistics as the size of the volume over which we 
calculate the star-formation density becomes small compared with 
the inter-particle distances. 

From Table [3] we can conclude that blow-out of gas in the 
models is not very efficient. Although the least massive models 
easily lose up to 80 per cent of their gas, this declines to around 
5 per cent for the most massive C models and for the D models. 
This difference in dynamical evolution is shown in Fig.[T2] show- 
ing the evolution of the gaseous component of 2 models, as a func- 
tion of time. While the C06 model is able to retain most of the 
gas and steadily converts it into stars, the C02 model is clearly 
not massive enough to be able to capture a large amount of gas. 
iMac Low & Ferr ara (1999) find that blow-away occurs for M gas up 
to 10 5 — 10 6 Mq , and that no blow-out of gas occurs for gas masses 
w 10 s — 10 10 Mq (depending on the energy released by SNe). 
These observations are broadly consistent with our results: gas loss 
is very small for initial gas masses > 4 x 10 8 Mq. Our upper limit 
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Figure 10. SFR (Mq /Gyr) as a function of time and distance from the 
galaxy centre, for four different models. The SFR is integrated over spheri- 
cal shells with width r max /50, with r max from the top to the bottom panel: 
3,2,1.2,0.6 kpc respectively. 
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Figure 11. SF density (MQ/Gyr/kpc 3 ) as a function of time and distance 
from the galaxy centre, for four different models. The spherical bins are 
equal to those of Fig. fTol 



12 S. Valcke et al. 




■1 le-2 4 



Figure 12. Evolution of the gas mass as a function of time (Gyr), 
for Models C06 (upper panel) and C02 (lower panel). The differ- 
ent lines show the fraction (mass divided by the initial gas mass) 
of the gas mass that is contained in a sphere with a given radius 
(lflo, 2Rc, 3R C , 4fl c , 5R C , WRe, 30 kpc). The top line shows the frac- 
tional gas mass remaining (i.e. not converted into stars.) R c is the half-light 
radius at the end of the simulation (see Table[5J. The top line for panel b) is 
not visible because it is ss 1. 



for blo w-away is fa 4 x 10 7 Mq (C01 model). Ferrara & Tolstoy 
2000) find a comparable limit for blow-out as Mac Low & Ferrara 



( 1999): Afgas ~ 10 — 10 Mq, also in good agreement with our 
results. 



4.4 Half-light radius 

The half-light radius, the radius of a galaxy containing half its to- 
tal luminosity, is an important diagnostic parameter in simulations 
of galaxy formation. Large deviations from values dictated by ob- 
servations hint at an erroneous spatial distribution of the star for- 
mation of model galaxies. Instead of using the three dimensional 
half-mass radius for R c , we chose to calculate the value accessi- 
ble to observers: the two dimensional half-luminosity radius. We 
first projected the light distribution of the models onto the xy plane 
and derived the cumulative luminosity function L(R) as a func- 
tion of radius. The value of R where L(R) reaches half its maxi- 
mum is then the half-light radius. Results for the models are shown 
in Fig. [13] where the logarithm of R a is shown as a function of 
Mb- Throughout the simulations we found that the half-light ra- 
dius of a model galaxy is very sensitive to the setup (e.g. initial 
conditions) of the simulations. Restraining it to a reasonable value 
is a far from trivial task, as can also be seen by the large varia- 
tions in the values of the CC02 models (Fig. 1 1 31) - Comparison of 
the data-points for the B and C models show that increasing the 
size of the DM halo results in an increase of the galaxy half-light 
radius. When the models are to be compared to the CC02 mod- 
els it should be to the CC02-B models, because those have initial 
densities closely resembling ours (the over-density of matter in the 
CC02-B models is set to 1: p gas = M gBS /Mx>M X Pu ~ 0.11 p u , 
whereas the CC02-A models have p g 



22 p u . We have p g 



M gas /M D M x 5.55 x 0.3 p u « 0.35 p u (see eq. (T5J)). Although 
somewhat too large, the model galaxies (B and C) perform reason- 
ably well in the range Mb € [—6, —15] mag. When we look at 
the larger models (A/b < — 15 mag) we see that R c decreases with 
increasing mass, signalling a much too central concentration of star 
formation compared with observations. A plausible explanation for 
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Figure 13. Logarithm of the half-light radius R c as a function of the ab- 
solute magnitude in the B-band Mb for several data-sets (observations: 
MA98 (Mateo 1998) (filled boxes: data th at only appears in this plo t), DR05 
JDe Riicke et alj 120051) . theory: CC02 IChiosi & Carrarol |2002|) . KA01 
(Kawata 2001)) as well as the models. 



this failure is that the breakdown of the models signals where the 
effects of hierarchical merging can no longer be neglected in the 
formation of galaxies, and thus that the simple initial conditions 
are no longer accurate in this region. When looking at the D mod- 
els, where we chose a for the DM haloes to be an arbitrary larger 
value (TableO, we see that it is possible to produce massive dwarf 
galaxies that do comply with observations, using a more extended 
halo than dictated by equation J22t . Here the correspondence with 
observations is excellent, showing not only that we can indeed re- 
gard the extent of the halo (within our set framework of star for- 
mation, feedback, ...) as the decisive factor concerning R c but also 
that we could fine-tune it, without changing any other parameters, 
to get excellent agreement between theory and observations over 
the entire dwarf magnitude range Mb £ [—8,-18] mag for R a . 
As the discussion in section l4~2l shows, however, these models de- 
viate from the expected Sersic profile. This deviation could be al- 
leviated by changing (reducing) a. We feel however that fiddling 
with parameters from model to model to be able to cover the en- 
tire magnitude range of the dwarf galaxies with excellent precision 
would not give much insight in the physics involved. A preferred 
strategy would be to study the evolution of dwarf galaxies in a fully 
cosmological setting, or to derive information about the initial con- 
ditions (extent and shape of the DM haloes, the ratio of gas mass to 
DM mass, ...) directly from numerical simulations. This is beyond 
the scope of this paper. We limit ourselves in this paper to the mass 
range where our chos en set of parameters give s good results. We 
note that the SAMs of INagashima et alj ( 120050 seem to be suffer- 
ing from the same behaviour of R c , it being rather large for small 
galaxies. There it is argued that this is a result of either the limited 
resolution of their N-body simulations, used to construct a merger 
tree, or that the SF timescale at high redshift is too long. 



4.5 Colour 

When comparing the B — V colour of the models with observations 
(Fig. I14t we find good agreement. Besides the value of B — V at 
t — 10 Gyr we also plotted the minimum and maximum value in 
the interval 7.4-10 Gyr for each simulation, as a function of Mb at 
those times (open tilted squares). This gives an idea of the variation 
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Figure 14. B — V colour of observed galaxies as well as the model galaxies 
as a function of blue magnitude M-g. . 
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Figure 15. Logarithm of the 1 -dimensional central velocity dispersion a 
(km s-1) as a function of the absolute magnitude in the B-band Mb for 
several data-sets, as well as the models. 



of B — V during the evolution of the model galaxies. This varia- 
tion is significant, which can be readily understood when looking 
at e.g. Fig.|7] as star-bursts tend to lower B — V because of the for- 
matio n of bright blue stars (see the discussion in lChiosi & Carrard 
2002). In this view the higher values for the CC02 models can be 
explained by looking at their fig. 2: all models have 1 predominant 
SF peak, no SF takes place after 6 Gyr of evolution, thus excluding 
the presence of high-luminosity young blue stars. Higher values of 
B — V could be obtained for the models by e.g. removing the gas 
content of a model after 1 or 2 SF peaks (as would be expected to 
happen in a cluster environment). When taking the state of Model 
C04 at 4 Gyr, which is between the first and second SF peak, and 
calculating luminosities as if the galaxy has an age of 10 Gyr, thus 
mimicking quiet evolution after the first SF peak, we get a B — V 
value of 0.71, close to the CC02 values. 



4.6 Central velocity dispersion 

We used the one dimensional central velocity dispersion, which is 
the three dimensional velocity dispersion (a) divided by \/3. Al- 
though, from an observational point of view, there are better mea- 
sures for a (e.g. light weighted average a) because of possible nu- 
cleation of dwarf galaxies, (Jid.c is a good enough measure for our 
comparison. We measured the central a by fitting an exponential 
function to the dispersion profile, of which we retain the central 
value. Results are shown in Fig. [15] Except for the least and the 
most massive C models, a is too high. The correspondence for the 
D models is however excellent. Overall we find reasonable conver- 
gence between models and observations within a factor of 2 (0.3 
dex). We further note that a comparison of Figs. 1 131 and 1 151 learns 
that raising a (going from the C to the B models) results in an in- 
crease of the half-light radius as well as a decrease of central a. 
Because of this we cannot increase a for the light models to get a 
better fit as this would cause R c to increase. The combination of 
R c and the central o thus provide an excellent tool to evaluate the 
soundness of halo properties within a certain framework of galaxy 
formation. 



4.7 Metallicity 

In Fig. [16] we plot metallicity (\og 10 (Z/ZQ)-B, with Zq = 
0.02), for the C and D models as well as 3 sets of data-points. 
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Figure 16. Logarithm of the metallicity of mode l galaxies (plotted at t = 
2.5,5, 1 Gyr) as well as data- points (NA04 JNagashima & Yoshiil2004l) . 
MIC07 iMichielsen et al. 2007)) as a function of blue magnitude. For the 
CC02 models the maximum metallicity is plotted. 



Z for our models is Zb, the light- weighted blue band metal- 
licity. New and improved information concerning the age and 
metallicity of dwarf galaxie s has become available during the re- 
cent years JGeha et al |[2003l: Ivan Zee et alj|2004l : IMichielsen et all 
2007; Pennv & Conselice 2007). Dwarf galaxies appear to be 
young, with [Z/B] approximately between -0.5 and 0, and a so- 
lar [a/Fe]. Three sets of points are plotted for each model, respec- 
tively at t — 2.5,5 and 10 Gyr. For the CC02 models the maximum 
metallicity is plotted. The overall trend is that metallicity increases 
with galaxy (star) mass (see e.g. Tremonti et al. 2004). We see from 
Fig. [T6] that the C models reproduce the observed trend well. The 
values agree reasonably with the data-points, being on average 0.2 
dex too high. For the D models the correspondence with the obser- 
vations is good. At 2.5 and 5 Gyr all the models correspond excel- 
lent with the data-points. 

In Fig. [T7] we plot Z as a function of stellar mass at 4 dif- 
ferent times for the C models. The initial decline of Z with stellar 
mass is a consequence of the rise of the number of SNe, which is 
not compensated for by the deepening of the galaxy potential well. 
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Figure 17. Logarithm of the metallicity Zb of model galaxies as a func- 
tion of the logarithm of star mass at 4 different times (from bottom to top: 
2.5,5,7.5,10 Gyr) for the C models. 



This is confirmed by Table[3] where the remaining gas mass is seen 
to drop when going from model C01 to C03. From the C03 model 
upwards, the metallicity rises with increasing galaxy mass, indicat- 
ing that the potential well is steep enough to compensate for the 
increase in energy input into the ISM. No clear trend concerning 
the remaining gas of these models can be seen in Table[3] 

This effect is not visible in the observational data. It is however 
unlikely that the decrease of metallicity with stellar mass would be 
observable in dwarf spheroidal galaxies as it depends strongly on 
the star formation in the range [3-10] Gyr (fig. [7}, and we expect 
that external processes (e.g. stripping) will greatly affect the evo- 
lution of these light systems. Furthermore this effect is model de- 
pendent and might be attenuated or even erased when modelling 
galaxies using a different DM halo. 



4.8 The fundamental plane 

The fundamental plane (FP) in K-space (Bender et al. 1992), as de- 
fined by the Virgo galaxies, is given by dBurstein et alj|l997l) : 



«3 



0.15 Ki +0.43. 



(24) 



where 



K\ = 



7f lo gio 



(Re 



n-2 



1 . J. 

= 7i loSlo( 



Re 



(25) 



K3 = 73 lo Sio(7rfc) 
I a is the effective surface brightness or the surface brightness 
within a distance R c from the galaxy centre. The physical mean- 
ing of the k variables can be described as follows: ki oc 
log 10 (M),K 2 oc log 10 (J e 3 A//L),K 3 oc log 10 (M/X). Our mod- 
els are plotted in «>space in Fig.[T8] We show an edge-on projec- 
tion of the FP in panel a), a face-on projection in panel b). In Fig. 
|19|we plot the FP using observed quantities (R c , a c , 7 C ) as coor- 
dinates. Fig. [20] shows the deviation of galaxies from the FP. The 
observational data on panel a) of Fig.[l8]learns that the trend of de- 
creasing M/L with decreasing mass is reversed at k,\ ~ 2.5. This 
indicates that dwarf galaxies and subsequent lighter systems suf- 
fer from depletion of baryonic matter relative to dark matter. Our 
models clearly reproduce this trend. From panel b) we learn that, 
despite a higher mass-to-light ratio, ki is lower for DGs and light- 
weight galaxies than for more massive galaxies on the FP. These 
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Figure 18. The fundamental plane in K-space fai, 1*2,03.). a) edge-on. 
The solid line is the equation of the FP derived from Virgo galaxies: 
k 3 = 0.15 ki + 0.43 lBursteinet"ai]|l997t corrected to Ho = 70 
km s —1 Mpc -1 ). b ) face-o n. Data represented by open squares are taken 
from Burstein et alj 119971) . data represented by asterisks are taken from 
iGeha et al.H2003l) . 



light-weight galaxies are therefore much fainter and more diffuse 
than their FP counterparts. Again we clearly reproduce this trend 
with our models. Figs.[l9]and|20]confirm this picture : our models 
reproduce the deviation of dSphs and dEs from the FP quite well. 
This deviation can hence be entirely contributed to the shallowness 
of the galaxy potential well, no stripping effects need to be invoked. 



5 DISCUSSION AND CONCLUSIONS 

We have presented N-body/SPH models of dwarf galaxies. The ba- 
sic SPH scheme was modified to include star formation, supernova 
feedback, chemical enrichment and gas cooling. The simulations 
start from a homogeneous gas cloud collapsing onto a DM halo. 
A first series of models (A models) was used to investigate the in- 
fluence of the gas particle number on the overall evolution of the 
system. We found that a minimum of 25 000 particles is required 
to resolve the evolution accurately. Using less particles leads to a 
reduced first star formation peak. Because of this, the ISM is in- 
jected with less energy (less SNe), resulting in higher subsequent 
SF peaks and hence a larger final stellar mass (Table [T). Based 
on these experiments, we decided to use 30 000 gas particles and 
15 000 DM particles in the simulations. 

All our model galaxies exhibit burst-like star formation. This 
is a direct consequence of star formation in a shallow potential well. 
Supernovae heat the gas, causing it to expand. The gas density con- 
sequently drops below the threshold p c for star formation. After the 
lapse of an adiabatic cooling period, the gas is again allowed to cool 
and contract. Star formation thus appears to be a self-regulating 
mechanism. Since substantial amounts of gas remain in and around 
the model galaxies throughout the duration of the simulations (see 
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Figure 19. The fundamental plane in physical coordinates: log 10 (.Re ) = 
-0.629 - 0.845 log 10 (J e ) + 1.38 log 10 (o- c ). The full black line shows a 
side projection of the FP. 
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Figure 20. Plot of the deviation from the fundamental plane. The black line 
represents zero deviation from the FP. 



fig- H2t . one expects that tidal interactions and ram-pressure strip- 
ping in dense environments will have a significant effect on their 
SFHs. 

We compared the model galaxies with observational data 
in terms of morphology (Sersic n, R c , no), kinematics (ctid.o). 
colour, metallicity and location with respect to the FP (k\ , K2 , K3). 
Despite their simplicity, the models do very well in the range 
—8 > A/b > — 15 mag. The models reproduce the slopes of the 
size, colour, metallicity, and velocity dispersion vs. luminosity re- 
lations very well. The zero-points of these relations are reproduced 
to within a factor of 2. The position of the model dwarf galaxies 



in the three-dimensional space spanned by (R e , ctid.o, ^e) agrees 
excellently with the observations. Dwarf galaxies trace a sequence 
within the FP (i.e. in the («i, K2) projection) that is almost orthog- 
onal to that of massive galaxies. In the dwarf regime, the height 
above the FP (in the (ki , K3) projection) is a decreasing function of 
galaxy mass. This is a consequence of less massive galaxies having 
shallower potentials than more massive ones, decreasing the ten- 
dency of the gas to fall into the potential well and increasing the 
capability of supernovae to support the gas by heating it. 

Furthermore we found that a good check for the validity of 
models is the combination of R c and <r c . R c is very sensitive to 
the location of star formation and hence to the time evolution of 
the density of the ISM. a c reflects, through the virial theorem, the 
total amount and distribution of (dark) matter. We demonstrate this 
using our B and C models, where it can be clearly seen (figs. [731 
and|15t that using a more diffuse DM halo indeed raises R c , and at 
the same time reduces ov. 
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